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The aim of this paper is to show how to extract dynamical behavior and ergodic properties from 
deterministic chaos with the assistance of exact invariant measures. On the one hand, we provide an 
approach to deal with the inverse problem of finding nonlinear interval maps from a given invariant 
measure. Then, we show how to identify ergodic properties by means of transitions along the phase 
space via exact measures. On the other hand, we discuss quantitatively how infinite measures imply 
maps having subexponential Lyapunov instability (weakly chaotic), as opposed to finite measure 
ergodic maps, which are fully chaotic. In addition, we provide general solutions of maps for which 
infinite invariant measures are exactly known throughout the interval (a demand from this field). 

Finally, we give a simple proof that infinite measure implies universal Mittag-Leffler statistics of 
observables, rather than narrow distributions typically observed in finite measure ergodic maps. 
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I. INTRODUCTION 

In many investigations we observe information that is 
generated by a physical system without prior knowledge 
of its nature. If such a system exhibits random-like be¬ 
havior, then we will be interested in constructing the 
probability density function from its distribution of data. 
It is well-known that even one-dimensional deterministic 
dynamical systems in the form Xt+\ = T{xt) can exhibit 
random-like dynamical behavior. In the case of so-called 
measuring-preserving maps, the density of trajectories at 
time t, i.e., pt.(x), has a peculiarity as t —> oo: if we it¬ 
erate randomly chosen initial points, the iterates will be 
distributed according to a limit density function p(x) al¬ 
most surely. In other words, pt(x) -A p{x ), where p{x) is 
the so-called invariant density of T, which satisfies 


/ p{x)dx = / p(x)dx 
Ja 


( 1 ) 


for any measurable subset A of the phase space. The so¬ 
lution of Eq. Q is the fixed point solution of the Perron- 
Frobenius operator pQ 


p(x) 


X 

4 i=T-\x) 


P(Zi) 

me,or 


( 2 ) 


where T' is the derivative of T and the sum extends over 
all preimages T^ 1 of the point x at which the density is 
to be evaluated. 

The concept of invariant density has implications that 
go beyond theory of one-dimensional maps. We can men¬ 
tion applications in random number generators, dynam¬ 
ical models for oil drilling and for the Hipp pendulum 
regulator, Poncelet’s closure theorem (a classical result 
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of projective geometry), the Bogoyavlenskij-Novikov cos- 
mological model (see Ref. [T] for a review of all these 
problems), and ways to characterize patterns of activity 
in the olfactory bulb [2], among other possibilities. As 
pointed out in Ref. [3], a further motivation for studying 
one-dimensional maps is that many higher-dimensional 
systems in the limit of strong dissipation approximately 
reduce to one-dimensional dynamics (Ref. [3] provides a 
variety of experimental evidences). 

Extracting statistical properties from deterministic 
systems based on long trajectories as data sets may not 
be a reliable procedure. It is possible that orbits display 
laminar or transient behavior for lengthy periods of time 
before settling into a fully chaotic mode. On the other 
hand, obtaining analytically an invariant density for a 
given map is typically a difficult task, when not unfeasi¬ 
ble. The procedure to get a deterministic map equation 
from an invariant density given beforehand has proved 
more fruitful in the literature, and is known as the inverse 
Perron-Frobenius problem (IPFP) [5H10j. For numerical 
computations of the invariant density, the most used and 
best understood algorithm is Ulam’s method mm, a 
procedure for approximating invariant densities by using 
matrix approximations of the Perron-Frobenius operator 
[Eq. § is the unity eigenvalue case]. 

Here we develop a new approach to address the IPFP 
based on the manner whereby monotone branches of a 
map are related to each other, an idea recently proposed 
in Ref. [13] to deal with the conjugacy problem of maps. 
Beyond providing exact solutions of the IPFP, the aim 
of this paper is to outline the main ergodic differences 
between maps of finite and infinite (non-normalizable) 
invariant measures. The invariant measure p is such that 

dp{x) = p{x)dx , (3) 

i.e., it is assumed as absolutely continuous with respect 
to the Lebesgue measure El. and will be the leitmotif 
for the results to be discussed here. 

Note that p, when infinite, is not a probability measure. 
For maps having infinite invariant measure the conven- 
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tional Lyapunov exponent vanishes, resulting in weakly 
chaotic behavior. Such kinds of maps are the motivation 
of so-called infinite ergodic theory |T5], an issue that is 
gaining increasing interest in the literature [HU - BIT] . De¬ 
spite its sound development, this theory still suffers from 
a paucity of models whose invariant measures are exactly 
known, and this paper aims to bring further alternatives. 
It is also noteworthy that, as far as I know, numerical 
methods for estimating infinite invariant measures do not 
seem to be available or hinted at in the literature E 30] . 
which highlights the importance of exact results. 

Exact invariant measures are also employed here to 
study the evolution of time averages of observables (like 
finite-time Lyapunov exponents) in a more efficient way. 
The approach, initially developed in Ref. [22], is ex¬ 
tended to more general shapes of maps, also enabling 
us to identify geometric characteristics of them. Based 
on this approach we also provide a simple proof that in¬ 
finite measure implies universal Mittag-Leffler statistics 
of observables, also the so-called Aaronson-Darling-Kac 
theorem, a central result of infinite ergodic theory [15] . 


II. MEASURE AND CONJUGACY 


Let us consider a Markov map T : I —> I defined on 
a partition {To,..., Ik- 1 } of / so that all of k branches 
T\j j = Tj are monotone. Thus, we can relate each branch 
Tj with the absolutely continuous invariant measure of 
the map in the form 

T j =*~ 1 o/z, (4) 

where 4>j are monotone functions. Note that each branch 
Tj is also a function extending on the whole interval I 
and with an image beyond I. Let us also introduce the 
foldings hj , relating each branch of T with T 0 , i.e., Tj = 
hj o To. Thus we have 

hj = TJ 1 o $, (5) 


where we simply set <I>o = $ and /io(x) = x. The cor¬ 
responding invariant density satisfying Eq. i> is such 
that 


poT 0 1 sgn {h'j) poT 0 1 oh j 1 

m ° To 1 1 ^ h'j o hj 1 | r 0 oTj'o hj 1 1 ’ 


( 6 ) 


where h'j and hj 1 denote, respectively, the derivative and 
the inverse of hj (the same holds for To), and sgn stands 
for the sign function. Within such an approach, the in¬ 
variant measure can be directly calculated by integration 
of Eq. |6|, namely [13], 


k-l 

p = $ + y^sgn(fe'-) $ o hj 1 . 
j=i 


(7) 


The choice of $ and hj enables us to construct a map 
equation with the desired absolutely continuous invariant 


measure. We can illustrate the use of Eq. ([7]) by means 
of a well-known unimodal map, the tent map 

T(x) = 1 — \2x — 11, (8) 

on I = [0,1]. This map was also considered by Lorenz 
as an approximation for the cusp-shaped Poincare first 
return map in the Lorenz strange attractor [31] . Here we 
have hi(x) = hj(x) = 2 — x and 4>(x) = x/2, resulting 
in the uniform invariant density p[x) = 1. 

A further advantage that we can draw from such an 
approach is to obtain conjugate maps. Two measure¬ 
preserving maps S and T are thus called if there exists a 
diffeomorphism oj such that [3] 152] 

S' o ui = u> o T. (9) 

Conjugacy gathers different maps into a same class where 
all maps share the same dynamics from the topological 
viewpoint. More specifically, conjugate maps have the 
same Lyapunov exponent and their invariant measures 
are simply related via uj: if T has measure p, then 

PS = P o w -1 . (10) 


Based on the folding approach above, the transformation 
u> relating both maps is such that m 

hs^ix = ujohj, (11) 

= ^oT 1 , (12) 


together with Cj +1 = Cj for all j, where ej = sgn (h'jh' s j). 
Furthermore, to is an even function if ey = —1. For exam¬ 
ple, we can ask for a map S conjugate with tent map (|8], 
let us suppose that hs{x) = x. Thus, the transformation 
is such that w(x) = ui{2 — x), with ei = —1. The solu¬ 
tion ui(x) = sn 2 (Kx, k) fulfills such conditions, where sn 
is Jacobi’s elliptic function, K(k) is the complete elliptic 
integral of the first kind, and k is its elliptic modulus [33]. 
By considering the addition identity of the sn function, 
this choice gives us Schroder’s map [32] ; 35j 


S K (x) 


4x(l — x)(l — k 2 x) 
(1 — k 2 x 2 ) 2 


(13) 


with measure ps = w -1 . The well-known Ulam’s conju¬ 
gacy u> = sin 2 (7rx/2) is the n = 0 particular case of map 
(13), resulting in the Ulam-von Neumann logistic map 
So = 4x(l — x) [IT] . 

Just for the sake of simplicity, we will focus only on 
maps with two branches and, therefore, our examples 
are either unimodal or piecewise expanding maps (see 
Fig. 0- Special attention will be given to maps having 
infinite invariant measures. For finite measure cases we 
also refer the reader to Ref. Ha¬ 


ul. INFINITE MEASURE MAPS 

Let us first consider interval maps, from [0,1] to itself, 
with a single marginal fixed point at x = 0: 

T(x -A 0) = 0, T'(x -A 0) = 1. (14) 
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The most common form of such maps is 
x t +i = f(x t ) modi, 


(15) 


with / : [0,1] —>• [0,2] bijective and h\{x) = x — 1, like 
the well-known generalization of the Pomeau-Manneville 
(PM) map [35] 


f(x) = x(l + x 4“), a > 0. 


(16) 


This kind of map was first obtained from Poincare sec¬ 
tions related to the Lorenz attractor for a = i EHES]. 

The remarkable characteristic of such maps is the inter¬ 
mittent switching between long regular phases (so-called 
laminar) near the indifferent fixed points and short irreg¬ 
ular burst ones elsewhere [see Fig. [l] (b)]. The applica¬ 
tions of such maps have been the most diverse, helping to 
model anomalous diffusion [37H40] , structure of natural 
languages m, vortex dynamics in evolutive flows [42] . 
weather systems [431 144j . DNA strands of higher eucary¬ 
otes 145] , and neural avalanches in the brain of mammals 
[iSj, among others. 

One problem to be faced here is the paucity of mod¬ 
els whose invariant measures are exact, i.e., analytically 
known for the whole interval. The most general result 
that we have so far is the relationship between the in¬ 
variant density and map equation near the indifferent 
point x = 0: p(x) ~ l/[f [x] — 1], up to a multiplicative 
constant; see Ref. [28] . 

By considering the approach discussed in Sec. |TT] Eq. 
5) is such that T 0 = / = <f> -1 op, and thus we have 
from Eq. (J?]) 

p(x) = $(x) + $(x + 1), (17) 

/( x) = 4> _1 [$(x) + $(x + 1)]. (18) 

The choice of $ gives us both the map equation and its 
corresponding invariant measure. Maps with marginal 
fixed points at x = 0 are such that 14>'(ar —>• 0)| —>• oo. 
Thus, from Eq. (18) one has 


/(a 


1 + 


<*>'(!) 

$'(a;)_ 


0, 


(19) 


provided we set <E>(1) = 0. 

Let us first consider the following choice: 


4 — 1 /“ - 1 


a > 0, a / 1, 


4>(x) = l 1 - !/« ' 

lnx, a = 1. 

For positive a^l we have Thaler’s map 


fa(x) = x 


1 + 


1+ X 


— X 


( 21 ) 


whereas a = 1 gives us the corresponding PM map ([16]). 
Note, by means of Eq. (19), that the behavior of Thaler’s 


map near x = 0 is the same of the PM map (16). For all 
a > 0, the invariant density is 


p(x) = b[x 


= h\^r- l l a 


(1- 


N-i/aq 


( 22 ) 


where b is an undetermined constant for the infinite mea¬ 
sure regime 0 < a < 1, since there is no possible normal¬ 
ization. Interestingly, Eq. (20) is the 1/a-deformed (or 
Tsallis 0?]) logarithm, i.e., $(x) = h ^/ a x for all a > 0. 
Thus, both a = 1 PM and Thaler maps are part of the 
same expression: 

f a (x) = exp 1/a [ln 1/a x + \n 1/a (x + l)], a > 0, (23) 


with exact invariant density given by Eq. (22), where 
expjy a x is the 1/a-deformed exponential. 

The a = 0 case also gives us a nontrivial class of maps. 
Two other infinite measure maps are suggested in Ref. 
[28], and we can also calculate exactly their map equa¬ 
tions and corresponding invariant measures for the whole 
interval. One model is 

$( 2 ;) = exp(a; - ^) — e, f3 > 0, (24) 

yielding the map 

fp(x) = ln -1/ ^{exp(aO' 9 ) + exp[(x + l)^' 3 ] - e}, (25) 

whose behavior near x = 0 is the same of the map em¬ 
ployed in Ref. [40] to produce strong anomaly diffusion. 
The other model comes from the choice 


$ q } r(x ) = exp(ln ? x r ) — 1, 


= 1 - 1/a, (26) 


where the parameter q mediates the strong laminar map 

(25) as q —> 00 and Thaler’s map at q = 1. For the choice 

(26) we have 


fq,r ( X ) — ®q [1 + ^ q,r( X ) 4 *&q,r( x 4 !)]• (27) 


For both cases (25) and (27) the invariant measure is 


given by direct substitution of corresponding $ into Eq. 

©■ 

From the models discussed above, many other maps 
with exact invariant measures can be generated via topo¬ 
logical conjugacy if we just choose a diffeomorphism oj 
from [0,1] to itself. For example, oj = f / 2 for any / given 
by Eq. (181 are possible choices. If we choose oj = f 1 / 2 / 2 
from map d2lb , a new rational map U = ui o Tj/ 2 0 w” 1 
with 


(20) conjugate with Thaler’s map T// 2 follows: 


U(x) = 


I f - x ( 1 + 2x ) 

2 /i /2 (2x) 1 + 2x _ 4x 2’ 


1 


fi /2 0 h\{2a :) = — 


x(l — 2x) 


6x + 4x 2 


x G [0,1/2), 


*€ [ 1 / 2 , 1 ], 


(28) 

with measure pu( x ) = P 1 / 2 0 The map U above 

has the same behavior of a = 1/2 in Thaler’s map near 
x = 0. 
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A seemingly different case from the indifferent fixed 
point maps discussed here is Boole’s map [48] 

B( x) = x -, (29) 

x 

from M\ {0} to itself. It is easy to see that map B(x) has 
uniform invariant density throughout a real line, and thus 
infinite measure, if we look at Boole’s integral identity 
f^ oo g{x)dx = J_ oo g o B(x)dx, valid for any integrable 
function g [48] , Despite not having an indifferent fixed 
point, Boole’s map behaves like a a = 1/2 map. For 
example, V = w _1 o B o to with u> = 1/(1 — x) — \/x on 
(0,1) yields the rational map [45] 


V(x) 


x(l — x) 


, x € [0,1/2), 


1 - 2x 
1 — 3x + x 21 


(30) 


se [1/2,1], 


with indifferent fixed points x = 0 and x = 1, and mea¬ 
sure = g-B ° w = w. 


IV. INVARIANT MEASURE AS A 
MEASUREMENT OF INFORMATION 

The exact characterization of an invariant measure is 
known to be a key issue in nonlinear dynamics because 
it rules the occupation probabilities of the iterates over 
entire phase space. Nevertheless, the predictive potential 
of exact measures goes beyond the statistical description 
of deterministic systems. We intend to explore here an 
idea that was recently introduced in Ref. |25j: the precise 
knowledge of the measure also enables us to estimate 
information by observing the transition dynamics along 
the phase space. 

We often turn to symbolic dynamics to simplify analy¬ 
sis of a dynamical system, representing their trajectories 
by infinite length sequences using a finite number of sym¬ 
bols. Since we are dealing here with two-branches maps, 
let us consider the usual partition of the interval [0,1] 
into two cells, A 0 = [0,2] and A\ = ( 2 , 1 ]. According to 
this approach, any trajectory {x t } can be represented by 
a sequence of binary digits {st} such that s* corresponds 
to the cell where x t belongs, i.e., s t = k for x t £ Ak, see 
Fig. 0 Then, redundancies that may appear in {s t } are 
eliminated by considering the algorithmic information, 
which is defined as the length of the shortest possible 
program able to reconstruct the sequence {st} on a uni¬ 
versal Turing machine |16j . 

Within the symbolic dynamics approach, easier pro¬ 
cedures are possible if the invariant measure is exactly 
known. To better understand how the invariant measure 
can assist in eliminating redundancies and making the 
calculation of observables more efficient, let us consider 
the following observable [251 : 


(31) 




FIG. 1: (Color online) Cobweb diagrams for unimodal (a) 
and piecewise expanding (b) maps. Instead of considering 
long trajectories as data sets, we translate the dynamics to 
states Ao,i. The shape of maps leaves specific features in the 
number of jumps from Ao to A ±, which can be inferred from 
exact invariant measures. 


where a(x £ A{) = 1 and a(x £ A 0 ) = 0. Equation (31) 


works like a filter, returning 1 whenever a transition from 
Ao to Ai happens and 0 otherwise. Now, we can compare 
the convergence of time averages for different observables 
under the same trajectories by using the Stepanov-Hopf 
ratio ergodic theorem [HI] 


EtptfP^aQ] 

EtJoVlTHx)) 


f ddg, 
J (fidg, 1 


(32) 


holding almost everywhere as t —> 00 for non-negative 
observables as long as they are integrable over /i. The 
limiting ratio (32) can also be regarded as a consequence 


of Birkhoff’s ergodic theorem [5T] for the particular case 
of finite measure ergodic maps. On the other hand, in 
the case of infinite measure ergodic maps [52j, we have 
to deal with averages of the type El=o 
which does not converge to a simple constant. In fact, it 
strongly depends on the initial condition x (more details 
in Sec. [Vb. Thus, the usefulness of ratio (32) is evident: 


it gives us pointwise convergence without depending on 
the finiteness of measure. 


X(x) = cr\T{x)\ [1 - ff(a:)], 
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By means of filter (311 we leave aside the complete de¬ 


scription of trajectories considering only the jumps from 
A 0 to Ai. From ratio (32) we have 


^tf[T fe ( 2 )] ~ 7 N t (x), 


7 = 


k—0 


I ddp 
f Xdn ’ 


(33) 


as t — > oo, where N t (x) is the number of jumps from Aq 
to Tli, given by 


t -1 

N t (x) = Y J X[T k (x)}. (34) 

k -0 

Thus, the calculation of f \dp plays a central role here 
since it depends on the partition z. For unimodal maps 
we have (see remark 0 ) 



MPo z<x«, 

lA T o7 1 ( z ), t i 1 ( z )\, z > x *, 


up to a multiplicative constant. Evidently, <j){x —> 0) 
rules <jr l (t —> oo). The generalized finite-time Lya¬ 
punov exponent is therefore given by d = In \T'\, but 
it is much more simple to calculate the number of en¬ 
trances Nt(x) than the summation Y^k=o \T' ( x )\ if d is 
exactly known. Thus, we have the alternative and more 
efficient formula 

A t (a;) ~ 7 ^/,^ ; t-t oo. (41) 


This approach was successfully employed for Thaler’s 
map; see Ref. [29J. Of course, Eq. (41) can also be 
used for finite measure maps and unimodal maps via Eq. 
(35) with ((t) = t. 


V. INFINITE MEASURE IMPLIES 
MITTAG-LEFFLER STATISTICS OF 
OBSERVABLES 


where 2 * is the nonzero root of x* = T(x*), whereas for 
piecewise expanding maps we have 



MpO 2 < O 

z > c. 


(36) 


It is easy to see that, for both Eqs. (35) and (36), the 
quantity f \dp has concave behavior with global maxima 
at x = x* and z = c, respectively. Thus, the coefficient 7 
has a unique global minimum at these partitions, thereby 
maximizing the number of entrances N t . Interestingly, 
the way how N t is maximized defines whether the map 
is unimodal or expanding, because 2 * and c have quite 
distinct roles in the evolution of trajectories. 

One of the applications is the calculation of finite-time 
Lyapunov exponents A t (x), which rules the asymptotic 
separation of initially nearby trajectories according to 


|<fai| ~ |$ 2 0 |exp[A t ( 2 0 )C(i)]- (37) 

For finite measure maps, it is widely known that ([(f) ~ t. 
On the other hand, for infinite measure maps, the most 
general instability regime is given by the sublinear growth 
rate 


C (t) ~ l a (t)t a , 


0 < a < 1 , 


(38) 


where l a {t) is slowly varying at infinity and with partic¬ 
ular asymptotes Z 0 (t 00 ) = 00 and li(t 00 ) = 0 . 
Slowly varying means that linr^oo l a (ax)/l a (x) = 1 for 
any positive a. The direct relationship between map 
equation and sublinear growth rate (38) was recently 
(see also Ref. 


given in Ref. 
C(t) ~ T 


1 


l (t ) 


), namely, 

sin(7ra)/7ra, 0 < a < 1, 
1/lnt, a = 1, 


(39) 


where T' ~ 1 — 1 /x(f>' as 2 —> 0. By comparison with Eq. 
(19) we have 

- r /i,/ m/ > a;^0, (40) 


x[x/$'{x)\' ’ 


Given that summations Z)L=o d[T k (x)} strongly de¬ 
pend on the initial condition 2 for infinite measure sys¬ 
tems, we are interested in the statistics of such quanti¬ 
ties when 2 behaves like a random variable. Following 
the Darling-Kac approach [52], the strategy for finding 
the asymptotic statistics is based on the calculation of 
moments, 


In(t) = E ■ 


E^( Tfc ( 2)) 


.k =0 


(42) 


with E{ } denoting the expectation value. 

The limiting ratio (32) tells us that the asymptotic 


statistics of J ra (t) does not depend on the choice of non¬ 
negative observable d provided it is integrable over p 
and that the summation T fe ( 2 )] is suitable nor¬ 


malized. Furthermore, ratio (32) also suggests the ex¬ 


istence of an asymptotically separable density function 
p(x,t) = Pt(x) ~ rtp(x) ~ r(t)p(x) such that 

t -1 . t -1 „ 

h(t) = XE / d{x)Pk{x)dx d(x)p(x)dx 

k= 0 d k =0 ’ 


r{t)dt / ddp 


(43) 


as t -> 00 , where r(t) does not depend on the observ¬ 
able, so that f r{t)dt cancels itself out. For Lyapunov 
exponents (d = In |T'|) we know that r(t) oc ('(t) due to 
subexponential instability, therefore p(x,t) ~ C'(t)p{ x ) 
as t — > 00 , up to a multiplicative constant. This asymp¬ 
totic separable form also proves to be quantitatively con¬ 
sistent with relationships between infinite measure and 
subexponential instability; see Ref. [251 . 

If, on the one hand, the limiting ratio (32) tells us 
that the statistics of observables is universal, on the other 
hand, Darling-Kac approach requires the use of a Markov 
process with stationary transitions. Given that T is a 
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Markov map, the choice of filter (30) via standard par¬ 
tition of the phase space meets this criterion since most 
part of laminar trajectories are compressed as “state A 0 ”. 
By performing the Laplace transform t > s on Eq. 


(43) with i) = Xi the first moment is given by 


h(s) f p{x,s) 


C(s 


sC(i 


x(x)dx —> C, s —> 0, 


(44) 


being C oc / yd/i a positive constant. The calculation of 
the second moment is as follows: 

m~2lf J x( x i)x( x 2)y2(xi,X2,t)dxid,X2, (45) 


where 


p(x 1 ,t 1 )p(x 2 ,t 2 -t 1 )dtidt 2 , (46) 

resulting in s/ 2 (s)[s£(s)] 2 —► 2!C' 2 as s -> 0. Similar 
calculations lead to the immediate extension for all n > 1: 


L(S) 

[<(*)]' 


n!C m , 


0. 


Now, by applying Karamata’s Tauberian theorem 
Eqs. (38) and (47), one finally has 


(47) 
H to 


lim E ■ 

t—> OO 


EtoxC^O*)) 


cm 


r n fl + a) 

= »l rn ' • (48) 

1(1 + an) 


For a = 1 we have the degenerate case 


StoXC^Qc)) 

m 


c, 


(49) 


which is similar to Birklroff’s ergodic theorem for linear 
m |5U- For 0 < a < 1, the numbers on the right 
side of Eq. (48) are just the moments of the Mittag- 


Leffler probability distribution with unit expected value. 
For a = 0 the Mittag-Leffler distribution becomes the 
exponential distribution 1 — e~^ on £ > 0. For 0 < a < 1 
the corresponding probability density function m Q (£) is 


m a (0 = 


F 1 /“(l + a) 

a£ 1+1 / c 


9o 


F 1 /“(l +a) 

{:l/a 


(50) 


where g a stands for the one-sided Levy stable density, 
whose Laplace transform is g a (s) = exp(—s“) [55] . Since 
the Mittag-Leffler distribution is uniquely characterized 
by its moments (48), one has 


lim Prob • 

t—¥ OO 


cm 


< z 


= / rn a {£)d£,. 
Jo 


(Si) 

The constant C is given by C = k J xdg, where k is 
such that C is well defined even in the infinite measure 
regimes; see Ref. |2S!j for more details. The quantity 


C((t) is known as a return sequence in the infinite ergodic 
theory [T51 . 


The results above give the statistics of the number of 
jumps N t but, as we have previously mentioned, such re¬ 
sults are also valid for any non-negative observables that 
are integrable over /i; see Eq. (33). It is worth mention¬ 
ing that the statistics of N t is usually considered in the 
literature by means of renewal theory. Despite provid¬ 
ing the correct Mittag-Leffler statistics of jumps, renewal 
theory does not provide the real number of entrances into 
A 1 ; for a detailed discussion see Ref. TH . 

Another important point is that the Darling-Kac ap¬ 
proach also tells us that, inversely, Mittag-Leffler statis¬ 
tics implies sublinear growth rate (38); see Ref. [53] for 
more details. 


VI. CONCLUDING REMARKS 


The goal of this paper was to investigate for which 
kinds of mappings there exist invariant measures, how 
many of them, whether they are equivalent (from the 
topological viewpoint), and what are their ergodic prop¬ 
erties. Typically, these issues are discussed in the frame¬ 
work of finite measure systems, separately from the infi¬ 
nite ones, which are relatively less known. Here we took 
a different perspective: by using a new IPFP approach, 
invariant measures are placed as the starting point, and 
their connections with the kind of instability (weakly or 
fully chaotic) and corresponding ergodic regimes are then 
drawn. For one-dimensional ergodic maps, Mittag-Leffler 
statistics of observables is an unequivocal signature of 
subexponential Lyapunov instability (weak chaos), and 
this is a direct consequence of infinite invariant measures. 
After some parametric change, e.g., from a < 1 to a > 1, 
if an infinite invariant measure becomes a probability 
measure, then Lyapunov instability becomes exponential 
and Mittag-Leffler density becomes a Dirac <5 function, 
which are signatures of finite measure ergodic systems. 
Thus, the strength of invariant measure plays the role 
of order parameter between weak (subexponential) and 
fully developed chaos. 

Subexponential instability implies zero Lyapunov ex¬ 
ponent in the conventional sense, and its relationship 
with infinite measures, as outlined here, may have par¬ 
allels to other kinds of systems. We can mention, for 
instance, chaotic attractors that lose stability in some 
invariant subspace by varying a parameter. For such sys¬ 
tems, the so-called blowout bifurcation [5B] occurs when 
the most positive normal Lyapunov exponent of the at¬ 
tractor crosses zero at the point of loss of stability. In¬ 
terestingly, at least two models exhibit infinite measures 
similar to a = 1 PM maps when undergoing blowout bi¬ 
furcation: chaotic motion of two identical dissipatively 
coupled one-dimensional mappings m and the drift- 
diffusion model [55[. 
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